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Abstract 

Analytical expressions for the crack-face displacements of an embedded elliptic crack 
in infinite solid subjected to arbitrary tractions are obtained. The tractions on the crack- 
faces are assumed to be expressed in a polynomial form. These displacements expressions 
complete the exact solution of Vijayakumar and Atluri, and Nishioka and Atluri. For the 
special case of an embedded crack in an infinite solid subjected to uniform pressure loading, 
the present displacements agree with those by Green and Sneddon. The displacement 
equations derived in this paper were used with the finite-element alternating method 
(FEAM) for the analysis of a semi-elliptic surface crack in a finite solid subjected to remote 
tensile loading. The maximum opening displacements obtained with FEAM are compared 
to those with the finite-element method with singularity elements. The maximum crack 
opening displacements by the two methods showed good agreement. 

Introduction 

Damage tolerant analyses require accurate calculations of stress-intensity factors at 
crack tips in two-dimensional and around crack fronts in three-dimensional bodies. Ex- 
perience with several crack configurations have shown that these cracks tend to grow 
in nearly elliptic, semi-elliptic and quarter-elliptic shape. Therefore, considerable atten- 
tion has been devoted to analytical and experimental studies of these crack configurations. 
While considerable data exists on the stress-intensity factors for these crack configurations 
[1-4], very little information exists on crack-face displacements. Crack-face displacements 
are directly measured in an experiment, while the stress-intensity factors are deduced from 
the experimental observations (such as stress freezing). Therefore, a need exists for ac- 
curate analytical crack-face displacements for comparison with experimental observations 
for three-dimensional crack configurations. Also, crack-face displacements are needed to 
develop three-dimensional weight function methods [5,6]. 

Recent literature has shown that the three-dimensional finite-element alternating 
method (FEAM) is a very economical and an accurate method for obtaining stress- 
intensity factors for elliptical and part-elliptical crack configurations [7-10]. The FEAM 
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uses the exact solution of an embedded elliptical crack in an infinite solid subjected to 
arbitrary crack-face tractions obtained by Vijakumar and Atluri [11], and Nishioka and 
Atluri [7,8]. While these solutions are comprehensive, the crack-face displacements are 
not available in references 7,8 and 11. In order to use the FEAM to calculate crack-face 
displacements, the exact solution of references 7,8 and 11 needs to be extended to derive 
the crack-face displacements for each of the polynomial pressure loadings on the crack 
faces. 

The purpose this paper is obtain the analytical crack-face displacements for an 
embedded elliptic crack in an infinite elastic solid subjected to arbitrary polynomial 
normal and shear tractions. First embedded crack solutions in the literature are reviewed. 
Next the fundamental equations relevant to the embedded elliptic crack problem are 
summarized. Then the crack-face displacements for mode I loading are obtained. Next 
the crack-face displacements for the mixed-mode loadings are derived. The displacements 
for the special cases of uniform pressure on the crack faces are compared with those from 
Green and Sneddon [12] and Kassir and Sih [13]. The crack-face displacements equations 
presented in this paper are used in conjunction with the FEAM to calculate the maximum 
crack opening displacements (COD) of a semi-elliptical surface cracks in a finite thickness 
solid. These COD values are compared to those obtained with finite-element method using 
singularity elements. 
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Symbols 

semi-major and semi-minor axes of the embedded elliptic crack 

magnitudes of the polynomial pressure loading on the crack faces defined by 

Eq. (1) and Eq. (47) 

elliptic integral of the second kind 

potential functions 

potential function defined by Eq. (7) 

elliptic integral of the first kind 

degree of truncation of the polynomials 

displacement along x a direction 

Cartesian coordinates with the origin at center of the elliptic crack with x\ along 
the major axis and X 2 along the minor axis of the ellipse 
shear modulus of the isotropic material of the solid 
(af - a\)/a\ 

a H a \ 

Poisson’s ratio of the material 
ellipsoidal coordinates defined by Eq. (3) 
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<r a p stresses in the solid 


Embedded Crack Solutions 

The exact solutions for an embedded crack in an infinite solid are the fundamental 
solutions used in alternating methods. Therefore, this section briefly reviews the classical 
works in the literature on embedded crack solutions. Some of the results from these 
references are used to build the present solution for the crack-face displacements. 

Exact analytical solution for a penny-shaped crack subjected to uniform pressure 
loading on the crack surfaces was obtained by Sneddon [14]. Later, Green and Sneddon 
[12] obtained the exact solution for an embedded elliptic crack in a infinite solid subjected 
to uniform pressure loading. Kassir and Sih [13] obtained an exact solution for an elliptic 
crack subjected to shear loading on the crack faces that can be expressed by a simple 
trignometric functions. Shah and Kobayashi [15] derived close form solutions for an 
embedded elliptic crack subjected normal loading that can be described by polynomial 
functions up to the third degree in terms of the two Cartesian coordinates that describe 
the ellipse. Smith et al [16] obtained the solution for a penny-shaped crack subjected to a 
polynomial pressure distribution on the crack faces. Smith and Sorensen [17] obtained the 
solution to an embedded elliptical crack in an infinite solid subjected to shear tractions 
that can be described by polynomial functions up to the third degree in terms of the 
two Cartesian coordinates that describe the ellipse. Later, Vijayakumar and Atluri 
[11] and Nishioka and Atluri [7,8] obtained the exact solution for an embedded elliptic 
crack subjected to arbitrary normal and shear loadings on the crack surfaces. Note that 
the solutions in references 12 through 17 are special cases of this general solution. As 
mentioned previously, the purpose of this paper is to extend the solution in references 7,8 
and 11 and derive the crack face displacements. 

Fundamental Equations For Mode I Loadings 

For completeness this section presents the equations that are relevant to an embedded 
elliptical crack in an infinite solid subjected to arbitrary tractions on the crack faces. The 
notation of reference 7 is followed. Some of these equations can also be found in references 
11-17. 

Consider an embedded elliptic crack in an infinite solid as shown in Figure 1. The 
crack faces are subjected to prescribed normal tractions <733 which can be expressed in a 
polynomial form as 
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( 1 ) 


11 Mm 

JO) - V'V' V x 2m-2n+i 2n+j 

ff 33 ~ / , / . / , 3,m-r»,n X 1 x 2 

i=0 j=0 m=0 n=0 

where * and j specify the symmetries of the load with respect to the axes of the ellipse, 
4?m-n,n are constants that give the magnitude of the load, and M is the degree of 
truncation of the polynomial. 

The stress distribution everywhere in the solid and the stress-intensity factors all along 
the crack border due to the prescribed pressure distribution of Eq. 1 has been obtained 
in references 7, 8, and 11. Here the crack opening displacements corresponding to the 
arbitrary pressure distribution of Eq. 1 will be obtained. 

Ellipsoidal Coordinates 

The elliptic crack perimeter can be described by 

(xi/ai) 2 -b (x2 /a 2 ) 2 = 1 , a l > a 2 (2) 

The coordinates x t (* = 1,2, 3) of any point in the solid may be expressed in terms of 
the ellipsoidal coordinates, £ a ,(a = 1,2,3) in the form 

A (A - A) x l = (A + Cl) (A + £2) (0? + £3) 


A (a| - a?) x\ = (a^ + 6) (<*2 + £ 2 ) (A + £ 3 ) 


( 3 ) 


A A x\ - & & £3 


where 


00 > & > o> ti > -A > £1 > -A 


Note that the ellipsoidal coordinates are defined by the roots of the cubic equation 


T 2 T 2 x 2 

f \ * X 1 x 2 x 3 

V(S) = 1 n— 

(Zj -f- S &2 *4“ 8 8 


( 4 ) 


= P{s)/Q{s) 


where 
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P(s) = (s - fl)(s - &)(» - 6), 


and 


Q(s ) = s(s-f aj)(s + a\) (5) 

In the X3 = 0 plane (the plane of the crack), the inside of the ellipse is given by 
£3 = 0, outside of the ellipse is given by £2 — 0, and the elliptic crack front is described 

by 62 = 6 = 0- 


Potential Functions, Stresses, and Displacements 

Using Trefftz’s formulation [7,8,11] for mode I problems, a potential function / 3 is 
defined such that 


where 


/3 = EE Fkl 

k l 


( 6 ) 


Fkl = 


d k+l 
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ds 


(V 


dxf dx[ J ( 3 y/Qis) 

In equation (6), are undetermined constants and cu(s) and Q(s) in Eq. (7) are 

defined in Eqs. (4) and (5), respectively. The partial derivatives of the potential function 
/3 are needed to obtain the stresses and displacements. They are denoted as follows. 


/ 3,0 

= F kt>p 

( 8 ) 


k 1 


h,Pi 

= Y1Y1 c w t'u.p'i 

( 9 ) 


k l 


f 3, pis 

= EE C 3 ,fc,/ F kl, P'16 

k l 

( 10 ) 


where is the partial derivative of fs with respect to xp y {/3 = 1,2,3) and so on. 
The stresses and displacements in the solid are given by 
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a\i = 2 /i[/ 3 ; ii +2 i// 3,22 + *3/3,311] 

(722 = 2 fl [ / 3j 22 +2j// 3) h + * 3 / 3 , 322 ] 

(712 = 2 n [ (1 - 2v) /3,12 + * 3 / 3,312 ] 

CT 33 = 2/x[-/3,33 + *3/3,333] 

( 7 3 i = 2^2:3/3,313 


( 11 ) 


and 


^32 = 2 /i * 3 / 3,323 

«1 = (1 - 2 1 /) / 3 ,1 + *3 /3,31 

U 2 = (1 - 2l/) / 3)2 + *3 /3,32 


( 12 ) 


«3 = -2(1 - v) / 3,3 + *3 /3,33 

where p is the shear modulus and 1 / is the Poisson’s ratio of the material. 


Crack Boundary Conditions 

The boundary conditions on the crack plane (X 3 = 0 ) are in two regions: inside and 
outside the elliptic crack. Inside the elliptic crack (£3 = 0 ), the boundary conditions are 


a (0) - 
ff 33 “ 


11 Mm 

EEEE A 

1=0 j=0 m—0 n=0 


2 m— 2n-ft 2n-fj 


(0) = a (0) _ 


'31 


32 


= 0 


(13) 


for 


(*l/oi ) 2 + (*2/<l2) 2 < 1, *3 = 0. 

Note that the right hand side of Eq. (13) are the prescribed tractions on the crack faces. 
Outside the elliptic crack (£2 = 0), because of symmetry in this mode I problem 
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(14) 


for 


U3 = 0 




(zi/ai) 2 + (*2 M) 2 > 1, *3 = 0. 

From the stresses in Eq. (11), the shear stress boundary conditions are identically 
satisfied and the normal stress boundary conditions require that 

1 1 Af TO 

-2 e/3,33 = EEEE us) 

t=0 y=0 m=0 n=0 

Because the potential function fz is expressed in terms of the unknown constants Cz t k,l 
and an unknown function F*/ (see Eq. (6)), Eq. (15) provides the relationship between 
(?Z,k,l and A? . Such a relationship was obtained from Eq. (42a) of reference 7. 

Equation (12) define the displacements everywhere in the solid. As the displacements 
at any point on the crack faces are sought, the displacements inside the ellipse (£3 = 0) on 
the 13 = 0 plane need to be evaluated. Because the derivatives Fjfcj^ (and hence /3,3-y) 
are finite when 13,^3 = 0 (see Appendix I for details), the crack-face displacements can 
be written as 


= (1 ~ 2i/) Jz,\ 


u 2 = ( 1 - 2 i /)/ 3i2 

U3 = -2(1 - v ) fz t 3 (16) 

Therefore, to obtain the crack-face displacements one needs the derivatives of the potential 
function fe with respect to xp. 


Mode I Crack Face Displacements 

Ffri t p must be evaluated to determine the partial derivatives of fe with respect to xp. 
The derivatives of F ^ are obtained by differentiating Eq. (7) as 


f m,p = 




dxf 1 dx Jts >/Q(s) 


roo | 


OJ 


«] 


k+l+l 


ds 


(17) 


where 
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k\ = k + Sip , /i = / + 8<ip , mi = 8zp 

and Sip, Sip, and 83 p are the well known Kronecker deltas. 

Because (^(£ 3 ) = 0, Eq. (17) can be written as 


F kl,p = 


00 


r °° d * 

J (3 dx\ x 




dxj dx™' 




Equation (19) needs to be evaluated in the interior of the ellipse (£3 = 0 ). Substituting 
oj(s) from Eq. (4) into Eq. (19) and performing the integrations one obtains 


k+l+l p q (_llP 

(fc + / + l)! ^2 / £ j J _ *1 h h *4 ts *6 Jp-q,q-r,r{( 3 ) (20) 


p=0 g—O r=0 


(2p - 2 g)! 

(P - ?)• 

(2g ~ 2r)! 
(g - r)\ 


2p-2q-ki 

U = {2p-2q-k l )\ 


h = 


2g-2r-/i 

X 2 

(2q-2r-h)\ 


*6 = 


2r-mi 



( 2 r-mi)! 


'6 (aj + s)P~<i [a\ 4- s)« r s r \A?( S ) 
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Equation (21) is an elliptical integral and can be written in terms of the Jacobian 
elliptic functions [7,18] as 


where 


Jp-q,<l-r,r — 


2 P +1 


f U (sn 2p f) (nd 2 *- 2r t) (nc 2r t) dt 

Jo 


r . r 


( 22 ) 


sn 2 u 


a l + & ’ 


dn 2 u = 


a 2 + fe 

°l + £3 ’ 


cn 2 u 


£3 

+ £3 


(23) 


The integral L P 9 _ rjr was integrated by parts and a recurrence relation is available 
in reference 18 (Eq. 357.01) and reference 7 (Eq. 30). A slightly different recurrence 
relationship that is useful for the crack-face displacement evaluation was obtained using 
integration by parts as shown in Appendix II. The recurrence relation is 


h.,-r,r = (2r [(»n»-‘u) (nc^—u) 

— (2p — 1) Lp_i ( g_ r —l t r“l “ (2? “ 2) Lp t q— rtr ^i] (24) 

Equations (19)-(24), provide all the information necessary to evaluate the partial derivative 
/ 3 ,/ 9 * 


Partial Derivatives With Respect To £3 

Because the crack opening displacements are sought the discussion is first focused on 
the derivative Fki y 3. To evaluate F^i 3, k\ = Jc, li = /, and m 1 = 1 need to be substituted 
in Eq. (20). 

On the £3 = 0 plane and outside the ellipse (£2 = 0), all terms in equations (23) and 
(24) are non-zero and the term x^ r ~ l in Eq. (20) is identically zero because r > 1 and 
£3 = 0. Hence, i^/ 3, and therefore U3 are identically equal to zero outside the ellipse, as 
required by the boundary condition of the problem. 

On the £3 = 0 plane and inside the ellipse (£3 = 0), because cn 2 u = 0 (see Eq. (23)), 
the integral L Vyq ^ ty r diverges. However, on this plane, £3 vanishes. The term 

2r — 1 r 
x 3 

takes an indefinite form and therefore needs to be examined when £3 = £3 = 0. 

The term is written as 
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*!'-* b.,-,r k £s =o = (2f _\) ^ H- 1 I- 1 "’-) 


- (2p - 1 ) z,' 1 ip-l,q-r-l,r-l “ ^ ( 2 9 2 ) x 3 ^' p,,-r ’ r ~ 1 ] i 3 =^ 3 =0 


= j [{J2-} 1 (sn 2p-1 u) (nd 2,-2r ~ l u) 

(2r — 1) k ' 2 l v cn u / 


-{2p-\)x\ r 1 L P - 1,^-r— l,r-l - « 2 (2g-2)x| r 1 L P.9-r,r-l] z3:= ^ 3:=0 ( 25 ) 

As will be obvious later, the last two terms in the square brackets of Eq. (25) are zero 
when X3 = £3 = 0. Now consider the term in braces in Eq. (25). With the help of Eq. (23) 

one obtains 


X S l 2r_1 _ T x 3 / 2 , c 

cnuJx 3 = 6=0 L VC 3 v 




T 2r— 1 

X3 = $3=0 


(26) 


— a 


2r — 1 
1 


^3 




2r— 1 

353 = ^ 3=0 


Noting that w(f 3 ) = 0 and utilizing equation (4), the term in the square brackets in 
Eq. (26) can be written as 

a\ a\ J 


Z3 


cn u 


2r-l 


= a 


2r— 1 


X 3 = ^ 3=0 


(27) 


Substituting Eq. (27) in Eq. (25) and using Eq. (23) gives 


2 1 4l 2r — 1 


, , 1 1 r_ f, x \ ?i\ 3 i 

c 3 L P,q-r,r I *3=6=0 ~ ( 2r _ 1 ) K i2q r 2 l 1 a\ a\ J J 


(28) 


Now consider the term x| r * Lp—i,q—r—l,r—l in Eq. (25). As before, considering the 
product of X3 and nc u terms, one has 


■>- 1 nc^-^u = ^[-S-l 

3 LcnuJ 


I 2r— 3 


(29) 


The term in Eq. (29) is identically zero when x 3 = f 3 = 0. A similar argument holds for 
the term x| r_1 L P)9 _ r>r _ i- 
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The partial derivative when X3 = £3 = 0 can now be written as 


F kl,Z 


113 


k+l +1 p q 

=e,=o = (k + , + 1)l £ (t+i+i-pV 


(- 1 )" 


<1 t 2 £3 t$ 


2 1 2r— 1 


1 2 1 1 [ xf X2HI 

(2 r - 1)! ■ a 2p+l ' (2r - 1) ' K '^ L 2 l a\ a\ i J 

Non-dimensionalizing Eq. (30) one obtains 


( 30 ) 


A+/+I p q 


F u, 31*3=6=0 - k+i rfi( k + l + 1 V- E EE 


(-1 ) p 


to as (* +/+1 - p > ! 


tl *2*3 


/ g l \2p— 2q— fe 2r— f 

1 <»1 7 ' l a 2 ' ' ( 2 p- 2 ?-Jfc)!( 2 g- 2 r-/)!( 2 r-l)!' 


1 

(2r ~ 1) 



*± 

«\ 


4\ 3r ^ 

a 2 


(31) 


Note that in the above equation the valid terms in the triple summation are those whose 
factorial arguments are either zero or positive integers. The opening displacements U3 are 
obtained as 


« 3 (*1,*2,0) = “2(1 -v) ^ C 3 ,Jfe,/.£fc/, 3 |x3=£ 3 =0 (32) 

k l 

The unknown constants are obtained in terms of the prescribed tractions on the 

crack faces by Eq.(42a) of reference 7. 

The non-dimensionalized opening displacements at the center of the crack (xi = %i = 
0) can be obtained from Eq.(32), where 


F kl , 3 I 


X 1 — X2 — X3 =0 


•HlJ+l 


(A: + / + 


( 2 * 2)1 ( 2 / 2 )! 
h I fc 2 ! / 2 ! 


Aa+Za+l ^jjAa+Za+r ! 2 r 
^ (fc 2 + h + 1 - r)! ' H* 2r — 1 

r— 1 

for even values of A: and / , with 2 k% = A:, 2/2 = /, and 


(33) 


F kl,3 \xi=X 2 — X 3 —' 0 0 


(34) 
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for odd values of k and /. Note that for odd values of k and l the antisymmetry of 

deformation requires that U3 be equal to zero along x\ and X2 axes, respectively. Equations 

(33) and (34) were evaluated for various values of k and / and Fkl,3 > s presented in Table 

1 . 

Note that Eq. (30) through (34) are equally valid for the penny-shaped crack (ai = 
a2 = a) because these equations do not contain any indefinite forms. 

Parial Derivatives With Respect To x\ and x<i 

The displacements ui and U2 (see Eq. (12)) require the evaluation of fz t \ and fz,2, 
and, hence, F^i i and F*/ 2 with *3 = £3 = 0. 

The partial derivative and are obtained by substituting 

k\ = k + 1 , l\ = / , mi =0 

ki = k , h = l + 1 , mi = 0 (35) 

in Eq. (20). This substitution yields terms like 

X3 2t L p>q -r,r (36) 

Evaluation of these terms can be performed in two parts, when r — 0 and when r > 1. 
First consider the later case. When r > 1, the term in Eq. (36) yields 

‘ ( 37 ) 

Because (x3/cnu) at X3, £3 = 0 is finite (see Eq. (27)), the term in Eq. (37) vanishes 
when X3 = 0. 

Therefore contributions to the displacements do not occur when r > 1. Now consider 
the term in Eq. (36) when r = 0. 

[4 r h,<l-r,r ] r=0 = W = -^+i J Q sn2p t nd 2q t dt (38) 

Instead of using the recurrence relation of Eq. (24), it is convenient to use ( see [18] and 
Appendix II) 

f sn 2p t nd* q t dt 

Jo 

where 
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I 2m = r nd 2m t dt (39a) 

Jo 

for 2m > 0 and 

I— 2m = G 2m = f U dn 2m t dt (39 6) 

Jo 

for 2m < 0. The integrals in Eq. (39a) and (39b) can be evaluated recursively by using 


T 2m(2 - k 2 ) l2 m + (1 - 2m) hm-2 ~ k 2 $nu cnu nd 2m+1 u 

W2 (2m + 1) ic* 5 


(40a) 


for 2m > 0 and 


^2m+2 


2m[2-K 2 )G 2m + (1 - 2m)/c ,S G 2m - 2 + k 2 sn u cn u dn 2m 

(2m + 1) 


(406) 


for 2m < 0. To evaluate the integrals in Eq. (39) and (40) one needs the starting values 
IoiGq, 1—2) and G 2 . These are obtained by setting m = 0 and m = 1 in Eqs. (39a) and 
(39b). This yields 


I 0 = G 0 = «, 1—2 = G 2 = E{ u) (41) 

As the displacements on the crack faces are sought, Eqs. (39) and (40) need be 
evaluated when 13 = £3 = 0. Thus, instead of incomplete elliptic integrals, one has 
complete elliptic integrals (with u = K(k)) and therefore Eqs. (40a) and (40b) reduce to 


for 2m > 0, 


hm+2 = 


2m(2 - k 2 ) I 2rn + (1 - 2m) hm-i 

(2m + 1) k ' 2 


G 2m+2 

for 2m < 0, and with 


2m(2 — k 2 ) G2m ± (1 - 2m)/ G 2m - 2 

(2m + 1) 


(42a) 


(426) 


lo = G 0 = K(k), I -2 ~ G 2 ~ E(k) (42c) 

where K(k) and E(k) are complete elliptic integrals of the first and second kind, 
respectively. 
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Therefore, the displacements u\ and U 2 at an arbitrary point on the crack faces are 
given by 


where 


«'y(*l.* 2 , 0 ) = ( 1 - 2 ^) ^2 C 3,k,l F kl,l\x 3 =U=0 

k l 


(43) 


F kln fU 3 ={ 3 =o 


*+/+i p 

(fc + / + 1)! E E 


p = 0 q — 0 


(~i) p 

(k + t+l-p)! 


(2p-2q)i (2q)l 

( P ~ ?)•' ?! 


2p-2q-ki 

X 1 


(2p — 2 q - k\)\ 


X 


2 q—li 
2 


(2g-/i)! 



(44) 


where 7 = 1 , 2 . Thus Eqs. (43) and (44) completely define the crack-face displacements 
ui and ti 2 . 

This completes the analysis of crack-face displacements for mode I type of deforma- 
tions. 


Mixed-Mode Crack Face Displacements 

For mixed-mode deformations, instead of a single potential function fa of Eq. ( 6 ) for 
the mode I case, there are three potential functions / a ,(a = 1,2,3) 

/« = EE F “ ( 45 ) 

k l 

where F^i is the function defined in Eq. (7). 

The displacements u a on the X 3 = 0 plane are given by [7,8,11-17] 

Ui (11,12,0) = —2(1 — v) /i ]3 + (1 - 2 i/) fz,\ 

« 2 (:ri,:r2,0) = —2(1 — u) /2 )3 + (1-2i/)/3,2 (46) 

«a(*l, 3 : 2 , 0 ) = -(1 - 2v) (fl t l +/ 2 , 2 ) - 2(1 — v) fz,z 
The first terms on the right hand side of Eq. (46) represent the displacements due to 
the mode II and mode III loading while the second terms represent the displacements 
due to mode I loading. 
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As in the mode I case, to evaluate the displacements one needs the partial derivatives 
fctfi (and hence Fkl t p) on the x$ = 0 plane and inside the ellipse (£3 = 0). However, all 
the partial derivatives are already obtained and are given by Eqs. (30) and (44). 

Thus for a general mixed-mode problem of an embedded elliptic crack in a infinite solid 
subjected to arbitrary tractions of the form 

11 Mm 

<’ 3 .(*l.* 2 . 0 ) = E E E E («) 

t=0 j —0 m=0 n=0 

where a = 1,2,3, the displacements on the crack surfaces are given by Eq. (46). The 
relationships between the magnitudes of A^ m _ n n and the constants C a / are already 
available from Eq.(42b) of reference 7 for mode II and mode III , and Eq.(42a) of reference 
7 for mode I deformations. 


Special Cases 


Elliptic Crack 

Consider an elliptic crack subjected to a uniform pressure of magnitude S on the crack 
faces. For this case k = l — 0 and A^’q 0 = S. Therefore, the opening displacements are 
computed using Eq. (32) as 


U3(*l . *2 , 0) = -2(1-1/) C 3 i 0,0 ^ 00,3 




(48) 


= -2(1-1/) C3 0,0 [ — — { 1 - - ~f } 

’ ’ loi 02 1 of a\ > 

The relationship between Aj^ 0 and C' 3^0 is [12,15] 


_ _2 ^ 3 , 0,0 
a 1 a 2 


2 /i 3,°,° 


or 


C-3,0,0 = Q1 j , S 


(49) 


(50) 


8 n E(k) 

Therefore, the opening displacements for an elliptic crack subjected to a uniform pressure 
S are obtained by substituting Eq. (50) into Eq. (48) and are 

«.(*!, «.o) = - (d/4) - (*!/«!)]“ (si) 
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which agrees identically with those due to Green and Sneddon [12]. The maximum opening 
displacement is 


u 3 (0,0,0) = 


(1 - v) Sa 2 


(52) 


/i E(k) 

Note that the total opening displacements between the crack faces is the COD and is 
equal to twice the value of the opening displacement U3(zi,-X2,0). 

Now consider the displacements tzi and «2- The displacements are 


— 4xi 


and 


2l/ ) ^3,0,0 *00 ,7 

(53) 

)0 >7 one obtains 


i i 

'IT 

i 

* 

i., i 

(54) 


(55) 


Equations (54) and (55) agree with those obtained by Kassir and Sih [13]. Thus the 
displacements u\ and 112 at any point on the crack face are 


ui( 


U2(xi,i2,0) = — (1 — 2u) 


5xi a 2 

\K(k) 

2/i af - 

Ie(k) 

5x 2 

1 

2/i a\ - a\ 

- K f2 


(56) 


goo 

E(k) 


Because of symmetry about the xi and x 2 axes, U2( I 1>0,0) = «i(0, X2,0) = 0- 
Eqs. (56) identically satisfy this requirement. The maximum and tt 2 displacements 
occur at the end of major and minor axes, respectively, and are given by Eq. (56) with 
z\ = d\ and X 2 — ct2* 


Penny-shaped Crack 

For a penny-shaped crack, ai = a2 = a, 2 £(k) — */2, an d Eq. (51) reduces to 

u 3 (xi,x 2 ,0) = 1 ~ ( x i /« 2 ) - fa!/ 02 )] * 

and the maximum opening displacement is 


(57) 
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(1 - v) 2 Sa 


« 3 ( 0 , 0 , 0 ) 


/* 


7T 


(58) 


Similarly, the displacements ui(n,i2>0) and U2( I 1 , 1 2 , 0) can be obtained from 
Eq. (56), using the limit as k — * 0 (see Appendix II) as 


ui(*l,*2,0) 


-(1 - 2 v ) 


Sx\ 

4/i 


«2 (^ 1 , 12 , 0 ) 


-(l-2i,) 


Sx 2 


(59) 


Use of the Crack-face Displacement Equations 

As previously mentioned, the crack-face displacement equations, Eqs. (30) and (44), 
can be directly used in the FEAM to calculate crack-face displacements for various three- 
dimensional crack configurations. To illustrate the use of these equations consider plate 
containing a semi-elliptical surface crack with ( 01 / 02 ) = 0.6. The plate is subjected to 
remote tensile loading as shown in Fig. 2. The crack configuration was analyzed by the 
FEAM and the three-dimensional finite element method with singularity elements using 
the models and methods described in referece 19. In the FEAM, the maximum degree 
of the polynomial ( M ) was equal to 5. The maximum crack opening displacement at 
( x\ — X 2 — X 3 — 0) obtained by both methods is compared in Table 2 for three crack sizes 
(defined by crack depth-to-plate thickness ratio) of 0.2, 0.5, and 0.8. The results from the 
two methods agreed well. For a deep crack, the differences between the two displacements 
is only about 3 percent. 


Concluding Remarks 

The evaluation of the crack-face displacements of an embedded elliptic crack in an 
infinite solid subjected to arbitrary tractions is addressed. The arbitrary tractions are 
assumed to be applied to the crack faces and are assumed to be expressed in a polynomial 
form in terms of the Cartesian coordinates that describe the ellipse. 

The exact solution obtained by Vijayakumar and Atluri, and Nishioka and Atluri was 
extended to obtain closed form expressions for the crack-face displacements. The crack- 
face displacements for special cases of uniform normal tractions agree identically with 
those obtained by Green and Sneddon, and Kassir and Sih. 

The evaluation of the crack-face displacements for embedded, surface, and corner 
cracked solids in finite bodies require the use of numerical methods such as finite element, 
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finite-element alternating and boundary-element alternating methods. The displacement 
expressions obtained in this paper were used with the finite-element alternating method 
(FEAM). An example of a semi-elliptic surface crack in a finite solid subjected to remote 
tensile loading was analyzed with the finite-element alternating method and finite-element 
method with singularity elements. The maximum crack opening dispalcements were 
calculated with both methods. In the FEAM, the crack-face displacements equations 
derived in this paper were used. The maximum opening displacements from both methods 
agreed well for both shallow and deep cracks. 
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Appendix I 


This appendix shows that the second partial derivatives, do not have singular- 

ities on the £3 = 0 plane, inside the elliptic crack (£3 = 0 ). 

The second partial derivatives F*/^ are [7,8,11,15] 


where 

TTlO 
with 

ki = k + l\ = / + ^2 7 ; mi = 1 + 8 (7.3) 

and 7 = 1,2,3. These derivatives are to be evaluated at £3 = £3 = 0 . 

The first term on the right hand side of Eq. (7.1) is given by Eq. (20) with k\, fi, and 
mi defined in Eq. (7.3). Using Eqs. (25) through (29), this term can be easily shown to 
be finite at £3 = £3 = 0 . 

Now consider the second term in Eq. (7.1). At £ 3 = £3 = 0 , this term can be written 
as 


Fki ^ - St 


00 ^l+Zi+mi 


k+l+l 


ds 


is dx ki dx l 2 ' dx? 1 ” y/Qffi 


+ F 


kl,3f 


( 1 . 1 ) 


M Ji _m 1 


(k + ' + 1)! ( 6 - 6 M 6 -W [ W“ W" % r ' ^ 

(1.2) 


<3, = (t + 1 + 1)! < 4 )l ' ' ( ¥ r ‘ In, fa-0 (M) 

Using Eq. (3) and rearranging, one obtains 


= (k+i+ix ££(=?)* ($'' ^ (fr-- 1 1-^=0 

Using Eq. (4), Eq. (7.5) can be written as 


1 


^3 pi - 1 




[1 - (xi /a x ) 2 - (x 2 /a 2 ) 2 ]s & " 


(/.5) 


m 
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For the terms 31 and F% 132 , m\ = 1 and hence = 1. Thus, no 

singularity exists in the right side of Eq. (J.6) and the term is finite. 

For the term F ^ 33 , mj = 2 and therefore Eq. (7.5) can be written as 


rpO 

33 


(* + / + !)! 


x l x 2 ( Zl \ k (—) 1 _ 1 _ 
a l ®2 a* l o| J V'S 


1*3, £3=0 


(1.7) 


This term has a singularity at £3 = 0. However, as shown in reference 15, when 
Eq. (1.7) is combined with the first term in Eq. (1. 1), the terms involving 1/y/lg cancel 
each other. Alternatively, because Fm is harmonic, 


^Jfe/,33 = ~^Jfc/,U ~ ^fc/,22 (^- 8 ) 

The derivatives 1 and ^fc|,22 can be eerily shown to be finite at 13 = £3 = 0. Thus 
the derivative F^i 33 also remains finite at 13, £3 = 0. 
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Appendix II 

The appendix derives the recurrence relationship of Eq. (23) and the integral in 
Eq, (39). This appendix also presents the identities of the Jacobian elliptic functions 
used elsewhere in the paper. 


Recurrence Relation Of Eq. (23) 


The integral L P)(? _ r|r is 

L„- r, r = [“ {in- 1 ’ t) (nd*"-' 1 ' t) (nc*' t) dt 

Jo 

Integration by parts was used to obtain the recurrence relationship as 
L Pjq - r ,r = I" {sn^-'t) {nd 2q ~ 2 t) {snt nc 2r t nd 2 ~ 2r t) dt 

Jo 

(nc2, " 1,) ( " <,2, " 2r " ,,) i: 

- (2p — 1 ) f (sn 2p ~ 2 t) {nd 2q ~ 2r ~ 2 t) {nc 2r ~ 2 t) dt 

Jo 

- k 2 (2q — 2) J“ ( sn 2p t ) {nd 2q ~ 2r t) (nc 2r ~ 2 1) df] 


Therefore 


(7/.1) 




h«-r,r = (2r Z 1 l) -72 M 2g ~ 2r - 1 «) 

— (2 p — 1) T'p-i.g-r- l,r-i — k 2 (2 q - 2) L Pi q- r>r -i j 

Derivation Of Eq. (39) 

Using one of the Jacobian identities the integral in Eq. (39) can be written as 


(J/.3) 
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f sn 2p t nd 2q t dt 

Jo 


L 


“ ('-f— F nd 2 ’ i dt 
0 « 


- JL Y' (~ t y P (p ~ ~ 2 ) ‘ ’ ( p i ± il f U dn 2j t nd 2q t dt 

_ «2 p 2^y '* j\ j 0 

J=0 


J_ y tlM. f u nd 2q ~ 2j 
K 2p Z_^ (p _ j)\j\ j 0 

j=0 


t dt 


1 y' (~0 J p ! 7 

_ k 2 p (p-j)!j! 2(9_j) 

In the above derivations and elsewhere in the paper the following identities of the 
Jacobian elliptic functions [18] were used. 


sn 


2 u + cn 2 u = 1 ; k 2 sn 2 u + dn 2 u — 1 


dn 2 u — K 2 cn 2 u = k/ 2 ; k/ 2 sn 2 u + cn 2 u — dn 2 u 
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cnu dnu ; d( nd u) 
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s/l — /c 2 sin 2 ^ 
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JO 
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Table 1 Values of at the center of the elliptic crack x\ = 12 = *3 = 0 
for various values of k and /, (0 < k t l < 5 ). 


k 

1 

a \a-lFkl,3 

0 

0 

-4 

1 

0 

0 

0 

1 

0 

2 

0 

32 /a? 

1 

1 

0 

0 

2 

32 /a! 

3 

0 

0 

2 

1 

0 

1 

2 

0 

0 

3 

0 

4 

0 

— 1536 /a? 

3 

1 

0 

2 

2 

- 512 /(a?a|) 

1 

3 

0 

0 

4 

— 1536 /a? 

5 

0 

0 

4 

1 

0 

3 

2 

0 

2 

3 

0 

1 

4 

0 

0 
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Table 2 Comparison of the maximum total crack opening displacement for 
a surface cracked plate by FEAM and finite element method with 
singularity elements. (ai/a 2 = 0.6) 


Crack Depth-to- 
Plate Thickness 
Ratio 

Mu 3 (0,0,0) E(k) 
(1-i/) 5 


FEAM 

FEM [19] 

0.2 

2.677 

2.699 


0.5 

3.132 

3.076 


0.8 

3.954 

3.839 
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